function realu = U1(x,t)
        s = 1/2*ones(size(x));
        ssign = 1;
        for k = 1:2:201
          s = s + 2*ssign*exp(-k^2*t)*cos(k*x)/k/pi;
          ssign = -ssign;
        end
        realu = s;
end
